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Abstract 

Massive colliding wind binary stars serve as labo- 
ratories for the study of strong-shock physics. In 
these systems fundamental flow parameters (veloc- 
ities, densities, directions) are often well known, 
and photon fields important for Compton scat- 
tering and ionization are usually defined or con- 
strained as well. In this paper we review observed 
X-ray and radio characteristics of some important 
colliding wind systems. We discuss how X-ray and 
radio observations of these systems can be used 
to test models of the collision of massive flows in 
the presence of intense UV radiation and magnetic 
fields. 

1 Introduction 



Though massive stars (stars with M > SMq) 
are relatively rare objects they have a cosmic 
impact which is disproportionate to their num- 
bers. They are chemical factories, they explode 
and they produce massive compact objects (neu- 
tron stars and black holes). Recent work sug- 
gests that the first stars which formed after the 
big bang had masses of tens or even hundreds of 
solar masses ( Abel. Anninos. Norman. &: Zhand . 
19981 ) and that massive stars may be the progen- 
itors of at least som e types of gamma-ray bursts 
( Heeer et al.l . 120031 ) . Over the last 2 decades the 
theoretical understanding of these stars has greatly 
advanced. However testing these theories is a chal- 
lenge since poorly-constrained quantities like con- 



vective core overshooting, mixing, weak-line opac- 
ity in the interior and ionization balances, opac- 
ities, instabilities and non-smooth density distri- 
butions ("clumping") in the atmosphere play key 
roles in determining stellar structure and evolu- 
tion. 

Historically observations of binary stars have 
played an important role in testing models of stel- 
lar evolution, interiors and atmospheres, and in 
testing Newtonian & non-Newtonian gravity as 
well. In this paper we expand this role by dis- 
cussing how binary stars can help test models of 
massive outflows and strong shocks. We discuss 
the symbiotic way in which radio and X-ray obser- 
vations of colliding wind binary systems help con- 
strain shock models, and how, by developing bet- 
ter understanding of colliding wind shock physics, 
we can develop a deeper understanding of massive 
star astrophysics. In ^ we discuss the properties 
of the stellar wind "laboratory", reviewing typi- 
cal properties of winds in different settings. In ^ 
we discuss the various types of interactions that 
can be usefully observed and modelled. We dis- 
cuss the fundamental physics of colliding stellar 
winds in ^ and show how they are applied to an 
important system, the extreme star rj Carinae, in 
33 and to WR 140, an important WR+0 binary, 
in ^ In 33 we consider how studies of collid- 
ing wind systems can help improve understanding 
of radiatively-driven mass loss. We discuss direct 
imaging of colliding wind shocks in 3H1 and how 
colliding wind emission can help define our under- 
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standing of the binary fraction in ^ Our conclu- 
sions and suggestions for future work are given in 

m 

2 Wind Characteristics of Massive 
Stars 

Stars with masses of > 8Mq possess suffi- 
cient radiative flux to produce dense mass out- 
flows from their photospheres by conversion of 
the outwardly-directed photon mome ntum to 
gas momentum by line abs o rption ( Mortonl . 
m^, ICa,stnr. Abbott, fc KleiuL Eli). These 
radiatively-driven stellar winds produce character- 
istic "P-Cygni"-type profiles in resonant transi- 
tions in the ultraviolet, in which outflowing ma- 
terial in front of the disk of the star produces 
blueshifted absorption, while material above and 
below the stellar disk produces predominantly red- 
shifted emission. Stellar winds are observationally 
characterized by three parameters: 

• the terminal velocity V^oo which represents the 
maximum outflow velocity of wind material; 

• the mass loss rate M which measures the di- 
vergence of material flowing out from the star 
and which is related to the mass density p by 
the conservation equation, M= 4:TTr'^V{r)p, 
where V is the wind velocity at any given ra- 
dius r; 

• and the wind acceleration parameter (3, which 
relates the velocity at a given radius to Voo as 
V{r)/Voo= (1 — R/r)^ where R is the stellar 
radius. 

In many instances and throughout this paper it 
is assumed that (3 = 1. The maximum blue ve- 
locity in saturated absorption components is an 
indication of the maximum speed of the outflow- 
ing stellar wind. Mass loss rates can be de- 
rived by detailed synthesis of the P-Cygni features 
since the density of the overlying wind determines 
both the optical depth of the P-Cygni absorption 
and also the strength of the emission component. 
However such fine analysis requires knowledge of 
the ionization balance through the wind, which 
is often poorly constrained. Mass loss rates can 
be determined more directly from observations of 



free-free emission in the radio, since the free-free 
brightness Si, in Jy is relate d to the mass-loss 
rate as 5^. oc {M /V^)'^/^ D''^ ()Wrig:ht fc Barlow . 
19751 ). where D is the distance to the source, 
assuming the wind is spherically symmetric and 
smooth (assumptions which are not well justi- 
fied). For a review of wind properties from a 
large selection of OB stars derived from UV and 
radio analyses, see Howarth &: Prinial ((l989|) and 
Bieging. Abbott, fc Churchwelll dlQsi i. 

Winds evolve as stars evolve. Typically an OB 
main sequence star has a mass loss rate in the 
range 10^'^ < M < 10^^ Mq yr^^ and wind speeds 
^ 1000 - 2000 km s~^ Theoreticallv Kx. ~ ^Vj^^c 
( Friend Sz Abbottl Il98d ) where Vesc is the escape 
speed at the stellar photosphere. As the star 
evolves the photosphere tends to expand, lowering 
the escape speed and thus lowering the observed 
wind speed, so that, in giant OB stars, the wind 
speeds are typically less than the winds of dwarf 
stars of similar spectral type. Because it's eas- 
ier to escape the photosphere of a giant star, the 
mass loss rates of giants are larger than those of 
dwarfs, typically by a factor of a few. Very mas- 
sive stars which enter the Wolf-Rayet phase exhibit 
very large mass loss rates (~ few xlO~^ Mq yr~^) 
with rather high terminal velocities. The most 
massive of these stars may enter the so-called "Lu- 
minous Blue Variable" (LBV) phase in which the 
star is thought to be near the Eddington stabil- 
ity limit. In this phase the star may experience 
episodes of extreme mass loss (~ few x 10~^ — 10~^ 
Mq yr~^) while the wind terminal velocities drop 
to a few hundred km s^^. 

3 Types of Interactions 



The wind from a massive star will eventually col- 
lide with something. A parcel of wind material 
can interact with another clump or stream in the 
wind. Stellar winds will collide with circumstellar 
material to form an evacuated region or bubble; 
this material may be primordial interstellar mate- 
rial or it might be material which was ejected from 
the star in an earlier evolutionary phase. In binary 
systems, the wind from one star will collide with 
the wind or surface of the companion. Some or 
all of these interactions can happen concurrently. 
In the following we briefly discuss the most com- 
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mon types of interactions which are known to give 
rise to observable X-ray and/or radio emission in 
massive stars. 

3.1 Self-Colliding Winds 

Radiatively-driven winds are inherently unstabl e 
to doppler perturbation ( Lucv &: Solomonl . Il97(l ) , 
since a small increase in the speed of a parcel of 
wind will move the parcel out of the shadow of the 
underlying wind material, producing an increase 
in the acceleration as the parcel sees more of the 
photospheric flux. This mechanism is thought to 
produce small clumps of shocked gas distributed 
throughout the winds of OB stars. These dis- 
tributed shocks are thought to be responsible 
for th e X-ray emission observed from single OB 
stars ( Berghoefer. Schmitt. Danner. &: Cassinelll 
19971 ). X-ray emission observed from OB stars is 
typically found to have T ~ 10^ K (suggesting that 
the shocks are not very strong) with little time 
dependence. Recent high-resolution X-ray spec- 
tra obtained by CHANDRA and XMM-Newton 
suggest that some of this emission arises about 2 
stellar radii ab ove the photosphere, or perhaps in- 



terior to this ( Cassinelli et al.l . l200ll : iKahn et alJ . 



20(1 ih . Such shocks may also produce observable 
non-thermal radio emission from single massive 
stars. 

3.2 Superbubbles 



Since massive stars evolve quickly, most do not 
have time to move very far from the nebulae out of 
which they formed. Most massive stars are found 
in HII regions and in clusters and OB associa- 
tions. These stars will directly interact with this 
rather dense environment by the transfer of wind 
kinetic energy into the circumstellar medium. This 
process can produce a wind-blown "super bubble" 
dCastor. McCrav. fc Weaver! . Il97.4 JDvsonl . Il977 : 
Weaver et al.l . Il9771 ) filled with hot shocked gas 
from the thermalization of the wind kinetic en- 
ergy. Such bubbles may be observable at X-ray 
and radio wavelengths. Since the wind proper- 
ties (wind speeds, mass-loss rates) vary as a mas- 
sive star evolves, the circumstellar region will in 
general be shaped by these various types of out- 
flows. For example, the fast thin wind of a main 
sequence star may carve out a large cavity within 



the primordial interstellar medium. If the star en- 
ters the LBV phase its massive, slow wind will 
flow into this cavity largely unimpeded. During 
the subsequent evolution of the star into the Wolf- 
Rayet phase, a less massive but faster wind devel- 
ops which can eventually overtake and collide with 
the slower wind from the earlier phase. All these 
collisions can give rise to observable X-ray and ra- 
dio emission, though important emission proper- 
ties (like abundances in the shocked material, the 
Mach number, the proportion of thermal to non- 
thermal emission, etc.) will generally be different 
at the different shock boundaries. 



3.3 Colliding Wind Binaries 

Perhaps the most interesting type of interaction 
occurs in massive binaries, where the wind from 
the primary star (here defined to mean the star 
with the stronger wind) can directly collide with 
the wind (or, in extreme cases, the photosphere) 
of the secondary star. This collision forms a "bow 
shock" around the secondary. Again thermaliza- 
tion of the kinetic energy of the wind means that 
portions of the bow shock can be heated to tem- 
peratures of millions of degrees and can radiate 
strongly at X-ray energies. This emission can be 
observable against the background emission from 
the "distributed shocks" in the stellar winds them- 
selves, since the colliding wind shock is expected to 
have temperatures about a factor of 10 or so higher 
than the distributed shock models if the wind- 
wind collision occurs where the winds have already 
reached their maximum velocities. In addition this 
emission should be highly variable due to intrin- 
sic changes in the emitting region, or to changes 
in the amount of intervening material between the 
observer and the shock, or to some combination. 
Along with X-ray emission these sources can pro- 
duce significant particle acceleration via accelera- 
tion of free electrons by the first order Fermi pro- 
cess. In the presence of a magnetic field these elec- 
trons will produce non-thermal synchrotron emis- 
sion in the radio. This population of accelerated 
particles could inverse-Compton scatter UV pho- 
tons to higher energies, producin g non-thermal X- 
rav o r even gamma-ray emission ()Eichler fc Usov . 
[l99.ih. 
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4 Simple Colliding Wind Physics 



Because the relative wind speeds are > 1000 km 
s~^, corresponding to a thermalization energy 
^ 1 keV, the shocked gas in a wind-wind 
colhsion will produce thermal X-radiation. 
The underlying physics which describes X- 
ray emission from wind-wind collisions was 
first discussed hv Ichexena.hchukI ^^ and 
by JPrihitskii fc Usovl h97(i) with ad ditiona l 



I 



199C) 



(1^), 
(1997), 



work bv iLuo. McCrav. fc Mac Low 
Stevens. Rlondiu. fc Pollockl (ll992lV lUso 

and mo r e rece ntly bv JPittard Sz Steven 

^^IHer] fl99^. and JHeulev. Stevens, fc Pittardl 
(|2003i). Though the shock structure can be 
complex and altered by numerous instabilities 
(see the paper by Julian Pittard at this meeting) 
the underlying physical model is actually fairly 
simple. Figur e U shows a cartoon of th e situation, 
adapted from IPrilutskii Sz Usovl ( 197a ). Consider 
a massive binary composed of 2 stars, star A 
and star B. In this system the dominant wind 
is determined by the value of the wind force 
ratio rj = pa/pb = {M aV^ a) / {M bVoo b) ■ If 
7/ > 1, the "bow shock" will be concave around 
star B; if it's less than 1, it will be concave 
around star A; and if ry = 1, the shock, in an 
idealized case in the absence of Coriolis forces, 
forms a plane midway between the two stars. 
The distance r^ from the center of star A to the 
bow shock, along the line of centers, is given by 
Tg = Dy/rj/{l + y/rj) where D is the separation 
between the two stars. If the shocked gas is 
adiabatic and the winds are at their terminal 
velocities, then the intrinsic X-ray luminosity of 
the shocked gas is L^; oc M /D. Phase dependent 
variations in the intrinsic X-ray luminosity will 
occur since the intrinsic luminosity depends on 
the separation between the stars, which varies 
with phase as long as the orbital eccentricity is 
7^ 0. The intrinsic X-ray emission from the shock 
is thus greatest at periastron and is at minimum 
at apastron. The observed luminosity is modified 
by any intervening absorption, L^^obs ot Lxe~''^" 
where Nh is the column depth to the X-ray 
source and a the appropriate weighted atomic 
cross-section. In general the absorbing material 
consists of interstellar medium absorption plus 
absorption in the wind of the frontmost star. The 




slow dejise wind 



hotter 
fast thio wind 



Figu re 1: Cartoon of a coll iding wind shock, adapted 
from IPrilutskii fc Usovl |l976i). 



component of wind absorption will generally vary 
as the stars revolve in their orbit. Thus, even 
when e = and the intrinsic X-ray luminosity is 
constant through the orbit, phase dependence in 
the observed emission can still be produced by 
phase-dependent changes in the wind absorption 
component for suitably high inclinations. 

Colliding wind systems also produce radio emis- 
sion by a variety of mechanisms. Particle accel- 
eration by first-order Fermi processes at the shock 
boundary can give rise to non-thermal synchrotron 
emission in the pr esence of magnetic fields near the 



shock boundary ( Doughertv et al.l . l2003h . Thus 



the presence of non-thermal radio emission may 
be an indicator of the presence of a wind-wind col- 
lision. Since the ionized stellar winds themselves 
are sources of thermal free- free emission, disentan- 
gling the contributions of thermal and non-thermal 
processes can be difficult. In addition observation 
of the non-thermal emission depends sensitively 
on the radio optical depth through the wind of 
the frontmost star, which is usually large for stars 
with large mass loss rates. In such systems the 
non-thermal component may well remain buried 
within the stellar winds, and unobservable. 

5 A Non-Canonical Example of Collid- 
ing Winds: r] Carinae 

The enigmatic massive star r/ Carinae shows emis- 
sion from various types of shock processes, from 
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ejecta-ISM interactions far from the star to vari- 
ations thought to be associated with a wind-wind 
colhsion near the star, rj Carinae is extremely lu- 
minous and believed to lie very near the Eddington 
Limit. The star has undergone at least two major 
eruptions since the 1840's. In addition to these 
ejection events (where the mass loss rate may have 
been 0.1 Mq yr~^ or larger) the star presently 
loses mass at a rate of about 10~^ — W^^Mq yr^^ 
and a velocity of Kx) ~ 500 km s^^(with perhaps 
higher velocities and larger mass loss rates near 
the stellar poles). The star is thus a laboratory 
to study the interactions of these ejections with 
the ISM, and the interactions of the wind with the 
ejections. 

The eruption of r/ Carinae in 1843 was associated 
with a large ejection of > 2Mq which now forms 
the bipolar "Homunculus" nebula around the star 
(see Figure [21). Although optical spectra recently 
found evidence of a high- velocity structure around 
the Homunculus which may be indicative of the 
forward shock just bevond the bipolar lobes 
( Dorland. Currie. Kaufer. &: Bacciottil . 120031 1 . no 
X-ray emission appears to be associated with this 
shock. Beyond the Homunculus lies the "outer 
debris field" identified as an incomplete optical 
shell which may be the remains of an e a rlier e rup- 



1971) 



or 



tion (|Walborn. Blanco. Sz Thackerav 
decelerated fast material associated with the 1843 
event. The outer debris field is associated with a 
"ring" of X-ray emitting material at the forward 
shock produced as the ejecta plows into the ISM 
(IChlebowski. Seward. Swank. Sz SzvmkowiaM 



IWeis. Corcoran. Bomans. Sz Davidsonl . 
2004f ). This outer shock is characterized by 



1984 



a temperature near 2 x 10^ K and absorption 
consistent with the ISM column as measured 
in the directions of neighboring stars. Figure 
m compares the X-ray and optical image of the 
nebulosities surrounding r] Carinae. 

5.1 Radio & X-ray Variations - Signs of a 
Companion Star? 

Recent evidence suggests strongly that rj Cari- 
nae is in reality a binary system with a 
5.53-year period. Radio and X-ray obser- 
vations in 1992 and 1993 showed what, at 
the time, was described as a "brightening" 




#. 



Ir»»^-H.pdkj 



Figure 2: Left: CHANDRA X-ray image of ry Carinae. 
This image was constructed from 5-100 ksec zeroth- 
order High Energy Transmission Grating Spectrometer 
observations. The colors represent different X-ray pho- 
ton energies. Red: 0.1 — 0.5 keV; green: 0.5 — 1.0 keV; 
blue 1.0 — 10.0 keV. The star appears as a white point 
source at the center of the image. Sharp linear fea- 
tures are either the dispersed spectra or instrumental 
readout strips. Right: HST/WFPC2 [Nil] 6583A im- 
age showing the nebulosity around t] Carinae; the star 
itself can't be seen in this image but is near the middle 
of the Homunculus. The scale and orientation of the 
optical and X-ray images are the same. 



or "outburst" of the star (IDuncari et al.l . Il995l : 
Corcoran. Rawlev. Swank, fc Petr j . |l995|)- This 
episode was associated with weakening of He I 
10830A emission and an abrupt change in the in- 
frared, ev ents which wher e sugg ested to be pe- 



and IPaminelil 



riodic by IWhitelock et al. 
()l99fil ). Continued monitoring at optical, radio 
and X-ray wavelengths showed that the emission 
over these wavebands is strongly correlated. Every 
5.5 years the radio, IR, optical and X-ray emission 
all experience a brief minimum in intensity. 

This multiwavelength variability has been inter- 
preted as arising from the presence of a hidden, 
hotter companion. UV flux from this companion 
is thought to play a significant role in modulat- 
ing the excitation and ionization of the nebular 
gas, which gives rise to variations seen in the He 
I 10830 line and in the free-free radio emission. 
By contrast the X-ray emission is believed to arise 
from the wind-wind collision and is thus thought 
to be a direct tracer of the conditions in the collid- 
ing wind shock. Figure El shows the observed X-ray 
brightness in the 2 — 10 keV band since 1996. The 
start of the X-ray minimum is given by 
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Time (Years) 




Figure 3: 2—10 keV X-ray brightness of r/ Carinae as 
measured by the RXTE satellite from 1996-2003. The 
X-ray "lightcurve" shows a gradual increase in X-ray 
brightness prior to the decline to a minimum which 
lasts 3 months. 



Column Densities 

• Chandra 

• XMH-Newton 

SAX (from previous cycle) 

• ASCA Ifrom previous cycle) 




Figure 4: Column densities derived from fitting the 
X-ray spectra obtained by ASCA, SAX, CHANDRA, 
& XMM-Newton. The column in the interval just after 
the start of X-ray eclipse is markedly higher than the 
column prior to the eclipse. 



UT X-ray minimum, years = 1997.95 -|- 5.536S 

(1) 

where the epoch is derived from daily monitoring 
observations with RXTE prior to the 1998 X-ray 
minimum, and the period, P = 5.526 years, is the 
interval between the start of the consecutive min- 
ima. The nominal uncertainty on the period is less 
than one day. 

Figure 0] shows the column densities, Nh, de- 
rived by fitting spatially-resolved X-ray spectra 
obtained by ASCA, SAX, CHANDRA, & XMM- 
Newton around the time of the X-ray minima in 
1998 & 2003. The column densities are clearly 
larger during the minimum than prior to the mini- 



Radio outburst of r] Carinae 




1992.0 1994.0 1996.0 



1998.0 2000.0 2002.0 2004.0 
Time 



Figure 5: The radio lightcurve of rj Carinae 
from Duncan & White from 1992-2004 (see 
[www . astro . umd . edu/$\sim$white/images/eta_time_f ull . html[ 



mum. In addition the interval of enhanced absorp- 
tion extends beyond the "recovery" of the X-ray 
flux. 

The radio lightcurve is shown in Figure El Though 
the frequency of monitoring of the radio emission 
is lower than that of the X-ray monitoring, a clear 
minimum in the radio flux separated by 5.53 years 
was observed. Unlike the X-ray emission the ra- 
dio emission reaches a maximum about halfway 
through the cycle, then begins a slow decline to 
minimum. Also, the radio flux just prior to the 
radio minimum apparently is declining with each 
advancing cycle, while the X-ray flux prior to the 
2003.5 minimum was substantially brighter than 
the flux just prior to the 1998.0 minimum. 

5.2 A Colliding Wind Model 

In the colliding wind model the periodic behav- 
ior of the X-ray emission is a consequence of or- 
bital eccentricity and the variation in the amount 
of absorption in front of the colliding wind shock. 
The X-ray emissivity and temperature distribution 
provide constraints on the mass loss rate and ter- 
minal velocity of the companion's wind, informa- 
tion which may be impossible to obtain in any 
other waveband. In turn, this information can 
be related to the photospheric escape speed using 
basic radiatively driven wind theory. This helps 
constrain the mass and radius of the companion. 
The photoionizing flux from the companion can 
be constrained from the free-free emission in the 
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radio. Thus, a nearly complete picture of the iden- 
tity of the hidden companion may be obtained by 
detailed modelling of the observed X-ray and radio 
emission. 

The X-ray spectrum of r] Carinae is charac- 
terized by a temperature of about 50 mil- 
lion K, which varies little throughout the cy- 
cle. This temperature corresponds to a pre- 
shock wind velocity of about 3000 km s~^, 
which is much higher than the measured 
wind speed of y Carinae itself, which is only 
500 km s~^ (JHilher. Davidson. Ishibashi. fc Gulll . 
200 ll ). so that this emission must arise in 
the shocked wind of the companion star. 



Corcoran. Ishibashi. Swank. &: Petrd ( 200lh com- 
pared simple colliding wind models to the RXTE 
X-ray fluxes. The orbit implied by this model- 
ing is shown in figure El This orbit is extremely 
eccentric (e ~ 0.9), with the companion behind 
7/ Carinae near periastron. Thus, in this model, 
the observed X-ray minimum corresponds to a 
time when the intrinsic X-ray flux is near maxi- 
mum intensity, but most of this emission is blocked 
from the observer by the wind from rj Carinae. 
The mass loss rates derived from fitting both the 
emission and the intervening wind absorption is 
Mri ~ 1^~^Mq yT~^ (which is about a factor of 
10 smaller than the mass loss rates attributed to 
r] Carinae based on radio and millimeter obser- 
vations) while the mass loss rate from the com- 
panion is Mc ~ 1O~^M0 yr~^. The wind charac- 
teristics of the companion (VJx),c ~ 3000 km s~^, 
Mc ~ IO^^Mq yr~^) are rather unusual compared 
to wind parameters derived from analysis of other 
m assive stars: in the studv of 24 Galactic O stars 
bv JRepolust. Puis. Sz Herrerd ( 2004 1 . the only star 
with similar wind parameters is rj Carinae's celes- 
tial neighbor HD 93129A (02If, Foo = 3200 km 
s-\ M = 2.6x1O"5M0 yr-i), though HD 93129A 
itself is a binary which might have biased the de- 
rived values. 

There are some obvious discrepancies between the 
observed X-ray emission and the colliding wind 
models: 

• for constant wind parameters the X-ray flux 
should be strictly periodic, yet the observed 
emission shows significant cycle-to-cycle vari- 
ations, which may indicate evolution of the 




Recovery 
Sep 23, 2003 

Recovery 

Mar e 1998 



10 20 

Distance (AU) 



Figure 6: Orbit of the companion star around 
77 Carinae based on the orbital elem ents from 
ICorcoran. Ishibashi. Swank, fc Petrel l)200l[) . 



mass loss from one or both stars; 

simple 2D models of the wind interac- 
tion suggest that the X-ray emission 
shoule be symmetric around periastron 
(IPittard. Stevens. Corcoran, fc Ish ibashi. 



19981 ). yet the X-ray brightness prior to 



eclipse ingress is about a factor of 3 higher 
than the brightness after recovery from 
the eclipse. This is perhaps indication of 
enhanced absorption caused by the distortion 
of the interaction region by coriolis forces 
through periastron passage (which will re- 
quire a full 3D model to determine). Alterna- 
tively. ICorcoran. Ishibashi. Swank, fc Petrj 
(|200l] ) suggested that this might be due 
to the enhancement of the mass loss rate 
from T] Carinae just after periastron passage, 
perhaps due to some interaction with the 
companion when the stars are close; 



unexpected variations 
on a timescale of 
Jlshibashi et al.', 'l997'; Ic 



or "spikes" occur 
80 - 100 days 



orcora n et al.l. 119971 : 
Davidson. Ishibashi. Sz CorcoranI 119981 ). 

The radio variation is more difficult to model be- 
cause the radio emission is dependent not only on 
the photon fiux from the companion, but also on 
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10 20 



Figure 7: O rbit of the WR star ar ound the 04 star 
adapted from lWhite fc Beckej l|l99,'il) . 



the distribution of the circumstellar gas, which is 
very non-uniform. The radio flux changes from 
an unresolved source during the X-ray/radio min- 
imum to a very extended source during the radio 
maximum, as significant regions of the Homuncu- 
lus nebula near 77 Carinae become ionized. The 
radio observations help constrain the amount of 
ionizing radiation from the companion star, which 
is key in determining the identity of the companion 
star. 

6 A (More) Canonical Example: 
WR140 



WR 140 (WC7-h04) is a long-period (P = 
2899 ± 1.3 days) , ecce ntric (e = 0.88) binary 
( Marchenko et al.l . 120031 ) . Like rj Carinae evidence 
of this system's duplicity was first obtained by the 
ide ntification of strong, p eriodic variations in the 



IR ( Williams et al.l . Il990^ ■ Unlike r] Carinae, the 



WR 140 spectrum shows absorptions which arise 
in the atmosphere of the 04 star so that radial ve- 
locity curve solutions exist for both the WR star 
and the companion. Thus the orbital elements and 
mass ratio are well constrained for this binary. The 
orbit of the WR star around the O star is shown 
in Figure [71 

6.1 X-ray and Radio Emission 

WR 140 shows non-thermal radio emission from 
the motion of shock-accelerated electrons in a mag- 
netic field near the wind-wind collision shock. A 
radio lightcurve for WR 140 thr ough its entire 8- 
yearcycle has been published bv lWhite fc Becker 
(|l995l l. They showed that the system exhibits 
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Figure 8: Radio observations of the 2. 6, fc 20 cm radio 
emission from WR 140, ad apted fromlWhite fc Beckeil 
f::i995l). including data from lWhite fc Becked l( 19 951) and 
IWilliams et alJ Ijl99ni ) , as a function of orbital phase. 
Phase = corresponds to periastron passage. 



strong variations at 2, 6, and 20 cm (Figure |H1), 
in that the radio emission at all wavelengths is 
greatest when the 04 companion is in front, and 
smallest when the 04 star is behind the WR star. 
This is presumably due to the the large opti- 
cal depth at radio wavelengths of the WR star 
wind, and th e smaller radio optical depth of the 
O star wind. IWhite &: Becker (1995 ) interpreted 
these wavelength-dependent variations to mean 
that most of the wind from the WR star is dis- 
tributed near the star's equatorial plane, noting 
that it was difficult to explain the observed maxi- 
mum at 20 cm when the O star is in front of the 
WR star, (orbital phase ^ « 0.85) since the wind 
from the O star should be opaque at 20 cm. 

WR 140 was recognized as an exceptionally bright 
stellar X-ray so urce in observ ations with the EX- 
OS AT satellite ( Polloclg, |l987|) . Subsequent obser- 
vations with ROSAT and ASCA and later RXTE 
showed that that X-ray emission also varied with 
orbital phase in a way roughly consistent with sim- 
ple colliding wind models: the observed emission 
reached a maximum near periastron passage, while 
the eclipse by the WR star wind of the X-ray emit- 
ting colliding wind shock shortly thereafter pro- 
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Figure 9: RXTE observations of tiie 2 - 10 keV X- 
ray em ission of WR 140, adapted from IPoUock. et alJ 
1|2004J) . Time in UT years is given along the bot- 
tom axis, orbital phase according to the ephemeris of 
iMarchenko et alJ 1)200,1) along the top. 




Orbital Phase 



Figure 10: Comparison of RXTE 2-10 keV X-ray 
emission variations of rj Carinae (black) and WR 140 
(in red). The rate of decline from the X-ray maxima to 
minima for WR 140 and rj Carinae, and the ratio of the 
maximum to the recovery level are similar in both stars. 
Interestingly a "flare" seems to occur just before the 
start of the decline to minimum in both stars. A data 
dropout is seen at « 2.09 in the WR 140 lightcurve. 



duces a brief minimum in the X-ray emission. The 
detailed X-ray variati on near periast r on p assage as 



observed by RXTE (|Pollock. et al.1 . 120041 shows 



• on p g 

Eoo^ 



a remarkable similarity to the X-ray eclipse of r] 
Carinae (compare similar phases in Figure |31 and 
in Figure ^ though in the case of WR 140 the 
eclipse is quite a bit narrower (see Figure ITUl The 
X-ray emission is thermal as revealed by ASCA 
and CHANDRA spectra though the hot emitting 
gas is probably o ut of thermal equilibr ium, at least 
near periastron ( Pollock, et al.l . l2004l ). 



7 Refining Winds 

Emission from wind-wind collisions depends upon 
wind properties so that determination of the level 
of emission can place severe constraints on mass 
loss rates and wind speeds. In addition emission 
from a wind collision zone can be used as a lo- 
calized probe of winds and can test how density 
varies with stellar radius, the smoothness of the 
stellar wind, and how the wind interacts with the 
radiation field. 

7.1 Clumping 

As discussed in §2, the most direct measure of the 
stellar mass loss rate is via radio observations of 
thermal free-free emission. Mass loss rates derived 
in this manner may be affected by the presence of 
non-thermal emission (produced, for example, in 
the wind-wind collision of unrecognized binaries) 
or if the wind is anisotropic or clumped. Clump- 
ing means that, rather than being composed of a 
smooth flow, the wind consists of dense regions im- 
mersed within a less dense outflow. Such clumps 
may arise due to the inherent doppler instability of 
radiatively-driven winds as discussed above. Mass 
loss rates calculated under the assumption that 
the wind is smooth and spherically symmetric will 
overestimate the true mass loss rate if the wind 
is clumped. Thus understanding how much of the 
wind is clumped is important to determining the 
actual amount of mass lost by the star, which plays 
a key role in determining how the star evolves. 

X-ray emission from the wind-wind collision in 
massive binaries can be used to help determine 
mass loss rates and the clumping fraction. Mass 
loss rates can be derived from the intrinsic lumi- 
nosity of the X-ray emission which is proportional 
to M^. The observed emission suffers absorption 
in the wind of the frontmost star, so that the ab- 
sorbing column to the X-ray source provides a 
measure of the density along the line of sight to 
the X-ray source. If the orbital geometry is known 
then the magnitude of the absorbing column to the 
X-ray source can be used to constrain the mass 
loss rate of the frontmost star by integrating the 
density profile through the wind of the frontmost 
star along the X-ray emitting portion of the wind- 
wind bow shock. Mass loss rate estimates from 
analysis of X-ray emission in colliding wind bina- 



10 X-Ray and Radio Connections http://www.aoc.nrao.edu/events/xraydio Santa Fe JVM, 3-6 February 2004 



ries has been done for 7 ^ Velorum ( Stevens et alJ . 
ll996':'Schild et all 120041) WR 14 (?), and rj Cari- 
nae (fPittard Sz Corcoranl . 120021 ) . In these cases 
the mass loss rates derived from the X-ray analy- 
ses are lower than those derived from estimates of 
the amount of fre e-free emission by a factor of a 
few. For example. IStevens et alJ (199a) derived a 
mass loss rate of 3 x 10~^ Mq yr~^ for the WC8 
star in 7^ Velorum based on an observation with 
the ASCA X-ray obsevatory when the companion 
was nearly in front of the WC8 star, a factor of 3 
lower than the mass loss rate derived from radio 
observations. More recently ISchild et al.l ( 20041 ) 
derived M« 2 x 10~^ from analysis of XMM- 
Newton observations when the O star was behind 
the WC8 primary (assuming the same distance to 
7^ Velorum, D = 450 pc, used by Stevens et al.). 
Schild et al.l ( 20041 ) showed that the difference de- 



rived from the X-ray absorption and the strengths 
of the X-ray lines imply that about 1/16 of the 
volume of the WC8 star's wind is clu mped. For 
T] Carinae, IPittard fc Corcoranl ()2002l ) derived a 
mass- loss rate from analysis of X-ray spectra which 
was lower by a factor of a few compared to the 
rates derived from other methods (like direct mea- 
sure of free-free emission and spectral synthesis of 
optical and UV lines), which might be indicative 
of either clumping, a non-spherical mass outflow, 
&:/or orbit-related variations in the mass loss rate. 



7.2 Radiative Interaction 



In colliding wind binaries the photospheric ra- 
diation fields can help mode rate the strength 



of th e wind-wind shock. IStevens &: Pollockl 
(|l994l ) considered the "radiative inhibition" of 
the stellar photons from the companion on 
the initial acceleration of the primarv's stellar 
wind , and vice versa. lOwocki fc GavlevI (|l995l ) 



and iGavlev. Qwocki. fc Cranmen (1993) consid- 
ered "sudden radiative braking" in which the 
winds can be almost completely decelerated by 
the UV flux from the opposing star. This might 
help explain why some short-period WR-star bina- 
ries like V444 Cygni have unexpectedly low X-ray 
fluxes. 



-1 1 r 




ThErmal WR147 wind 




Figure 11: 



WR147 imaged by MERLIN 
[ I1997I1 . The southern source is resolved 



/illiams et a,lJ 

thermal emission from the wind of the WR star, while 
the northern source is non-thermal emission formed in 
a coUiding-wind bow shock. The crosse s mark the po- 
sitions of the optical stars identified by iNiemela et al.l 
{199^. Right: CH ANDRA HRC X-ray image from 
IPittard et all (|2002|) . 



8 Direct Imaging of Colliding Wind 
Bow Shocks 



The ability to resolve sources of emission in col- 
liding wind systems is an important tool to help 
understand the physics of the emission mecha- 
nisms. A case in point is the star WR 147, 
a nearby Wolf-Rayet star which was resolved 
bv VLA dChurchwell etHI . Il992l ) and MERLIN 
( Williams et al.l . |l993) into a southern and north- 
ern source separated by 0.6" (Figure ITU) . The 
southern source was identified as thermal free-free 
emission originating in the wind of the WR star. 
Interestingly, at 5 GHz this thermal source is not 
spherically (or circularly) symmetric. The north- 
ern source shows non-thermal synchrotron emis- 
sion. The emission is produced by electrons which 
are Fermi accelerate d in a colliding wind shock. 
Niemela et al.l ( 1998h identified an 0-type star just 
behind this shock front, revealing WR 147 as a 
long-period colliding-wind binary. 

Pittard et al.l ( 200j) attempted to resolve the col- 
liding wind shock from WR 147 at X-ray ener- 
gies using the CHANDRA High Resolution Cam- 
era, which offers a spatial resolution of about 0.5". 
This observation revealed that the X-ray emission 
is clearly extended but the HRC lacked the spatial 
and spectral sensitivity needed to clearly distin- 
guish emission from the WR star (or O star) wind 
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from the bow shock. 

9 Finding Binaries 

It's important to understand the fraction of mas- 
sive stars which are binaries in order to understand 
how massive stars form and evolve. However con- 
ventional means of determining the duplicity of a 
star can be of limited use in determining whether 
or not a given massive star has a companion. Bi- 
naries can be identified by periodic radial velocity 
variations in their photospheric absorption lines, 
but in massive stars such lines are often rotation- 
ally broadened or contaminated by emission or ab- 
sorption in stellar winds. In systems with large 
mass ratios the more massive star will show only 
small-scale radial velocity variations. These prob- 
lems are especially acute in long-period systems 
which usually show little radial velocity or photo- 
metric variability for most of the orbit. 
In principle colliding wind emission can be use- 
ful for identifying massive binaries, especially 
long-period systems where the winds collide at 
terminal velocities and where t he sh ocks are 
strong. lOoughertv fc Williams! (|2000|) consid- 
ered using non-thermal radio emission as a 
tracer of the massive star binary population. 
DouEfhertv fc Williams! (EooO) examined 23 WR 



systems, found 9 non-thermal radio emitters, and 
showed that 7 of these 9 non-thermal emitters 
were known to have a spectroscopic or visual bi- 
nary companion. They also showed that not all 
binaries exhibited non-thermal emission, however. 
In particular short-period WR binaries tended to 
be thermal sources, presumably because the non- 
thermal emitting region is buried within the opti- 
cally thick WR wind. 

Cherepashchuk! ( 197a) suggested using X-ray emis- 
sion from wind-wind collisions to determine the 
binary fraction of WR stars. As is the case 
with non-thermal radio emission, stellar wind ab- 
sorption may play a decisive factor in determin- 
ing how much of the intrinsic emission is observ- 
able, especially at energies E < 2 keV where 
wind opacity is greatest. Sources which emit at 
E > 2 keV will suffer little absorption for even 
large column densities {Nh < 10^^ cm~^), so 
in this sense X-ray observations have an advan- 
tage over radio observations at detecting collid- 



ing wind emission. However the rather coarse spa- 
tial resolution of X-ray imaging versus the ability 
of radio imaging at milli-arcsecond levels is im- 
portant. Foi_exam2le_^ROSA21.ajid_j4^ 
vations (jRho. Corcoran. Chu. fc Reach! . !200l! l de- 
tected a bright source of hard X-rays associated 
with HD 164492A, the 07 star which ionizes the 
Trifid nebula, and suggested that this star might 
be a colliding wind system. Subsequent higher 
spatial reso lution observations with CHANDRA 
( Rhoet al.! . Enn.1) showed that HD 164492A is ac- 
tually a soft source and is not overluminous in X- 
rays compared to single stars, and that the hard 
emission detected by ROSAT and ASCA is from 
an unrelated nearby object. Because massive stars 
usually congregate in associations and clusters, it 
remains to be seen how many of the catalogued "X- 
ray overluminous" massive stars may be resolved 
into two or more sources in higher spatial resolu- 
tion observations. 

10 Conclusions 



Because colliding wind emission can be identified, 
and the parameters on which it depends (Voo, M) 
observationally quantified, observations of collid- 
ing wind binary systems at X-ray and radio wave- 
lengths allow for tests of stellar wind structure, 
and for understanding shock processes and modes 
of instability. Colliding winds in massive bina- 
ries provide good laboratories for testing models 
of shock-generated radio & X-ray emission. C om- 
parison of simple analytic models ( Usov! . !l992h and 
more complex hydrodynamic simulations (Pittard, 
these proceedings) show that gross characteristics 
of the emission can be reproduced. In general, ob- 
servations of both radio and X-ray observations 
are necessary to fully test models. Radio obser- 
vations provide high spatial resolution imaging to 
probe the spatial distribution of the shock struc- 
ture, and can also provide a measure of the mag- 
netic field strength near the shock. X-ray spectra 
can provide measures of the spatial distribution of 
absorbing material in front of the shock, and can 
be used to measure the flow dynamics by deter- 
mining the temperature distribution and also via 
high resolution studies of X-ray emission lines. 

New high resolution X-ray and radio observations 
however often show complicated behavior which is 
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not yet well matched by the models. This sug- 
gests the importance of realistically modelling the 
effects of spatial and temporal variations in the 
stellar winds and Coriolis forces in order to deter- 
mine the observed emission properties. The com- 
plex nature of the wind flows and their apparent 
inhomogeneities indicates the important need to 
develop realistic models in 2 &;/or 3-D and the im- 
portance of including magnetic fields. 
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